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Sound can be just as quantum as light. But our toolbox for single quanta of sound, i.e. phonons, 
is currently insufficient. Here we describe a new component that enables a chip-based, solid-state ana- 
logue of cavity-QED 1 2 utilizing acoustic phonons instead of photons, phonitons instead of polaritons. We 
show how long-lived and tunable acceptor (hole) impurity states in silicon nanomechanical cavities can 
play the role of a matter non-linearity for coherent phonons just as, for example, the Josephson qubit 
plays in circuit- QED 3 . This system enables the control of single phonons and phonon-phonon interac- 
tions, dispersive phonon readout of the acceptor qubit, and compatibility with other nano/optomechanical 
components 4 5 such as phonon-photon translators. Phonons, due to their unique properties, enable new 
opportunities for quantum devices and physics. 



Circuit-QED has revolutionized the field of cavity-QED 
(cQED) providing a stable platform for light-matter inter- 
action in the microwave regime along with large coupling 
strengths and solid state integrability. Progress in the field has 
enabled potential applications such as single microwave pho- 
ton source^ and quantum logic gates 3 on a chip. In an ideal 
crystal environment, phonons may play a role analogous to 
photons, though they propagate with the much slower speed 
of sound. That acoustic phonons can be quantum coherent 
has been explored in a number of novel architectures, allowing 
seminal experiments in nanomechanical cooling^ElE] trapping 
of phonons in artificial phononic bandgap cavities 5 9 , pho- 
ton translation via phonons 1011 , and indirect qubit-phonon 
coupling 12 13 . What is missing to complete the analogy for 
phonons is a non-linear element similar to an atom in cQED. 

Such an element has been proposed recently 14 where a sin- 
gle donor in silicon (P:Si or Li: Si) couples directly to confined 
phonons to form a hybrid state: a phoniton (direct analog of 
the polariton in cQED). The impurity -phonon interaction may 
be large due to large deformation potential matrix elements: 15 
(\\f s f\Dij\\\f s ) ~ eV. The proposed system, however, requires 
very high frequencies (a few hundred GHz) and can be diffi- 
cult to integrate with other phonon components. While other 
impurities such as in diamond^ or in III-V semiconductors 
are possible, a practical system in silicon would be highly de- 
sirable given recent demonstrations of high-Q cavities in sili- 
con nanostructures 9 , silicon's investment in materials quality, 
and compatibility with CMOS technology and silicon photon- 
ics. 

In this Letter we propose a new quantum circuit element 
based on a single acceptor (such as B, Al, In, etc.) embedded 
in a patterned silicon nano-membrane, that would be compat- 
ible with optomechanical components 5 9 . The acceptor qubit, 
Fig. 1, is easily tunable in the range of 1 — 50 GHz by external 
magnetic field, and also by additional electric field or strain. 
We show that the regime of strong resonant and also strong 
dispersive coupling of the qubit to a confined acoustic phonon 
is possible. We give explicit measurement protocols to ob- 
serve the phonon vacuum Rabi splitting and readout of the 
phonon number state fine structure. 

Engineering the qubit levels. Two acceptor qubit arrange- 
ments are possible based on the lifting of the 4-fold ground 
state degeneracy via external fields (see Methods). For a mag- 




FIG. 1. (a) Hole valence bands in Si; 4-fold degeneracy at the 
band top (and of lowest acceptor states) corresponds to particles of 
spin J = 3/2 (Tg representation of cubic symmetry, see, e.g. RefP^). 
(b) Ground state splitting via external magnetic field along [0,0, 1] 
direction; allowed (forbidden) phonon transitions and qubit phonon 
driving (text). Level rearrangement is via additional strain. Sys- 
tem manipulation via electric static/microwave fields is possible, (c) 
Examples of nanomechanical ID and 2D phonon bandgap cavities 
reminiscent of already fabricated high-Q cavities in a patterned Si 
membrane 5 9 . The acceptor is enclosed in the cavity and an on-chip 
phonon waveguide allows phonon coupling in/out of the system. 



netic field H z = (0,0, H z ) along the crystal [0,0, 1] growth di- 
rection one can choose the lowest two Zeeman levels, = 
1 3/2), \§2) = 1 1/^)' as the qubit, which is the primary fo- 
cus of this paper (Fig. lb). The Zeeman type interaction is 
given bjW H H = M X JH + g' 2 (/ x 3 ^ + c.p.)}; here c.p. 
is cyclic permutation of jc, y,z\ J x -> • • etc. are the spin 3/2 ma- 
trices (in the crystal directions), and the renormalized g-values 
g[, g' 2 (Mb is Bohr magneton), depending on the acceptor 
bound states, fulfill the relations \g[\ « 1, \g' 2 \ <C \g[ fBEEl . 
The qubit energy splitting 8£# ~ Vog[H is tunable in the 
range ^1—40 GHz for H = 0. 1 — 3 T. The qubic term in Hh 
lifts the level equidistancy: the lowest and highest splitting 
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(Fig. lb) are larger than the middle splitting by 



81 



0.09. 



For a field H tilted away from the crystal axis [0,0,1], the 
qubit splitting is weakly angle dependent. 

Alternatively, a second qubit arrangement is possible with 
mechanical stress in addition to the magnetic field. Mechan- 
ical stress lifts the ground state degeneracy only partially: 
e.g., for stress along the crystal growth ^-direction (Fig. lc), 
states I ± 3/2) and | ± 1/2) remain degenerate. Providing the 
stress causes an energy splitting larger than the splitting due 
to magnetic field H z , the level configuration in Fig. lb rear- 
ranges so that the lowest (qubit) levels will be 1^) = | — 1/2), 
\§ f 2 ) = 1 1/2). This forms an alternate qubit, decoupled from 
phonons to first order; however, coupling can be switched on 
via electric field (see below). The effect of strain is described 
by the Bir-Pikus HamiltonianEl 
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The renormalized deformation potential constants a 1 \b' \d' for 
acceptor B:Si can be estimated in the effective mass approx- 
imation (EMA) 15 , or extracted from experiment 17 ; the lat- 
ter gives bf = -1.42 ± 0.05 eV, d' = -3.7±0.4eV. Using 
H £ , we estimate a splitting of 5£" £ wl - 10 GHz for exter- 
nal stress in the range of ~ 10 5 — 10 6 Pa. Such stress can be 
created in tensioned nanomechanical structures/membranes, 
improving also the mechanical g-factoF^J A much larger 
stress (> 10 7 Pa) is possible due to a nearby (random) crys- 
tal defect 20 or for [0,0, 1] strain due to a SiGe substrate. The 
resulting splitting in the few hundred GHz range effectively 
suppresses the qubit-phonon coupling. 

Strong coupling of acceptor to confined acoustic phonon. 
In order to calculate the qubit-phonon coupling (Methods) we 
consider a quantized phonon field in addition to any classical 
field. The coupling of the qubit { 1 3 /2) , 1 1 /2) } , is calculated to 

a plane wave mode 8 vac ^ e~ iq r with polarization and 
energy hv a q, that proved to be a good estimation of coupling 
to modes with realistic boundary conditions 14 (Fig. lc). The 
relevant matrix element is proportional to the "phonon vac- 
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v c is sound velocity, and V is the mode volume). For long- 
wave acoustic phonons, the coupling is (the acceptor is placed 
at maximum strain of the nanomechanical cavity): 
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where the polar angles of the wave vector q are with respect 
to z-direction. Thus, the mode ti has a maximum along the 
phonon cavity (0 « 7C/2), Fig. lc. An alternative is to have 
an in-plane magnetic field H x along the crystal [1,0,0] in- 
direction (the latter is chosen to be along the phonon cavity): 
both modes t\ , ti (now at « 0) are preferably coupled to the 
cavity. The maximal coupling gmax,cy 2 scales as y/ayjj V 
yfqj V, as expected for a (Is — )> Is) transition. Taking a cavity 
volume V ~ dX 2 (d = 200 nm is the Si membrane thickness) 



we get coupling in the range g/2% ~ 0.4 — 21 MHz for the fre- 
quencies of 1 — 14 GHz (Table 1). The other allowed transi- 
tion 1 3/2) — » I — 1 /2) (at twice the qubit frequency) is well de- 
tuned, while the transitions 1 3 /2) -> | - 3/2) , 1 1 / 2) -> | - 1 / 2) 
are phonon forbidden (Fig. lb). 

Generally, when the in-plane magnetic field has some an- 
gle 0o with the cavity (crystal x-axis), all transitions are al- 
lowed. Also, the qubit coupling to a preferably confined cav- 
ity phonon will change. As a qualitative example, consider a 
plane wave transverse mode t\ (or ti) along the x-axis (0^0): 
the coupling will change in the same way as in Eq. ([2]), with 
replaced by 0o. This allows manipulation of the qubit-cavity 
coupling by rotation of the magnetic field. 

Acceptor qubit relaxation rate. The qubit relaxation in the 
cavity is bounded by the bulk phonon emission rate r re i ax < T, 
and reads: 

ri2(9o) = 2tefe K( cos2 2eo+1) [2/3vf + l/v, 5 
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2/v? + 3/v?]}; (3) 



here the contribution of longitudinal phonons is only a few 
percent (speeds of sound in Si: v/ ~ 1.7 x v t = 8.99 x 
10 3 m s _1 ) and results in Table 1 are for 0o = 0. Note that the 
coupling in this case can be switched off (e.g. for a ^-mode 
along x-direction, at 0o = 7c/2) while the relaxation cannot. 

For the alternate qubit, { | — 1 /2) , 1 1 /2) } , the stress and mag- 
netic field are parallel along the growth z-direction (Fig. lc). 
Here both coupling and relaxation are zero in the absence of 
electric field and can be switched on using non-zero elec- 
tric field E z in the same direction (Methods). It can be 
shown that the qubit-phonon coupling is given by the same 
Eq.Q multiplied by a coupling factor, a function of the 
splitting ratios r h = r e = ^: f(r h ,r e ) = (^z+Z- - 

l)/V(l+z+)(l+z_), with z± = (1 ± V(l TrW+r 2 T 
r h) 2 / r e- Thus, e.g. for = 0.5 — 0.9 this factor reaches 
w 0.25 — 0.65 for some optimal value of the electric field split- 
ting, r e <l, which allows strong qubit-cavity coupling. 

The relaxation time (Table 1) is comparable to T\,T2 
measurements 23 in bulk Si at low B:Si doping (8 x 10 12 cm -3 
or 500 nm acceptor spacing), where 7i ?exp ~ 7. 4 /us and 
r 2 * exp ~ 2.6 jus. While bulk T 2 * is limited by electric- 

dipole interacceptor coupling 24 , both T\, T2 may improve for 
defect- free samples^^, isotopically purified Si 25 , and single 
acceptors 20 ; T\ may also improve in nanomembranes (d <C X) 
due to phase-space suppression. Finally, the relaxation to pho- 
tons is strongly suppressed since (v t /c) 3 w 10 -15 . 

Phonon cavity loss. In the lD/2D-phononic bandgap Si 
nanomembranes considered here (Fig. lc), the main cavity 
loss mechanism is due to (fabrication) symmetry-breaking ef- 
fects, coupling the cavity mode to unconfined modes, and also 
due to cavity surface defects 26 . Bulk loss mechanisms are 
negligible in the GHz frequency rang e u \ 26 [ The cavity Q- 
factor for frequencies of one to few tens GHz can be engi- 
neered in the range of 10 4 — 10 5 , or even higher^a. 

Collecting all calculations in Table 1 , one finds that strong 
resonant coupling is possible: g\2 ^> Ti2,K in a wide fre- 
quency range, allowing a number of Rabi flops, ~ 30 — 100. 
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parameter 


symbol 


circuit-QED 


Quant Dot-QED 


B:Si(l GHz) 


B:Si (4 GHz) 


B:Si (8 GHz) 


B:Si(l Tesla) 


resonance freq. 


CO r /27T 


5.7GHz 


325 THz 


1GHz 


4 GHz 


8 GHz 


14GHz 


vac. Rabi freq. 


g/27T 


105 MHz 


13.4 GHz 


0.41MHz 


3.27 MHz 


9.26 MHz 


21.4 MHz 


cavity lifetime 


1/k,G 


0.64 /us, 10 4 


5.5ps, 1.2 10 4 


15.9^/s, 10^ 


4/us 


2 /is 


l.U/us 


qubit lifetime 


l/r 


84 ns 


27 ps 


386.5^/s 


6 /is 


0.75 /is 


0.14/us 


critical atom # 


2TK/g 2 


< 8.610" 5 


% 1-87 


<4.91(T 5 


<21(T 4 


<3.910" 4 


<6.910" 4 


crit. phonon # 


r 2 /2 g 2 


< 1.6 10" 4 


<9.410" 2 


<5.110" 7 


<3.21(T 5 


<2.610~ 4 


< 1.410- 3 


# Rabi flops 


2g/(K+r) 


-98 


-0.8 


-79 


-99 


— 64 


-34 


cavity volume 


V 






0.037 X 5 


0.148X 3 


0.296 X 5 


0.52 X 5 


wavelength 


X 


5.26 cm 


921 nm 


5400 nm 


1350nm 


675 nm 


385 nm 


dispersive coupling 


X = 8 2 /^ 


17 MHz 




0.04 MHz 


0.33 MHz 


0.93 MHz 


2.14MHz 


peaks' resolution 


2x/r 






- 199 


-25 


-9 


-4 


# of peaks 


2X/K 


-70 




-8 


- 16 


-23 


-31 



TABLE I. Key rates and parameters for circuit-QED 21 (ID cavity), Quantum dot(QD)-QED 22 (InAs QD embedded in a GaAs photonic 
crystal) vs. the B:Si phoniton system in a patterned Si membrane (of thickness d = 200 nm) phononic bandgap cavity; we show calculations 
for maximal coupling at frequencies of 1GHz, 4 GHz, 8 GHz, and 14 GHz, for cavity volume V — dX 2 and Q— 10 5 . The limiting frequency 
for strong dispersive coupling is reached at « 21 GHz, when % = T. Number of Rabi flops, dispersive coupling, and dressed state's resolution 
parameters, 2%/T, 2%/k, are comparable to that in circuit-QED and significantly surpass semiconductor quantum dots. 



Here, a frequency of 1GHz is the limit for T ~ 20 mK un- 
less an active cavity cooling is performed 9 . Different energy 
scaling of the coupling and qubit relaxation sets a maximal 

energy splitting: E™ x = ^ \J « 200 GHz. At such high 
frequencies the limiting factor would be the significant dete- 
rioration of the cavity g-factor. Still, e.g., at 14 GHz even 
Q = 10 3 leads to strong coupling. 

Observing the vacuum Rabi splitting. A suitable observable 
could be the averaged phonon cavity field amplitude \(b)\, 
which we have calculated, Fig. 2a-d, in the so-called "two- 
state approximation'^. On Fig. 2a we show the field ampli- 
tude spectrum for co fl = co r = 8 GHz as two Rabi peaks at ±g 
vs. external phonon driving frequency COj, while Fig. 2c shows 
the "anticrossing picture" of the spectra at different acceptor- 
resonator detuning A ar . Here the left resonance width and 
(normalized) height are Tl = sin 2 0or + cos 2 0oK and cos 0o, 
respectively, while for the right resonance they are obtained 
via sin 0o <-> cos0o, with tan(20o) = 2g/A ar . 

Another experimental option to observe the dressed reso- 
nance^) is to sweep the qubit detuning (via the H -field) while 
keeping the probe input in resonance with the phonon cavity, 
Fig. 2d. As seen from Fig. 2c, the field amplitude will exhibit 
a resonant dip; it is shown for four different cavity frequen- 
cies, co r = 1, 4, 8, 14 GHz. For large detuning the resonance 
is approximately a Lorentzian (Fig. 2d) with a full width at 
half maximum, {fwhm} « £ 2 /k (a weak dependence on the 
qubit relaxation F is suppressed for strong couplings). 

With increasing temperature the Rabi peak width will be 
broadened 28 by the factor 1 + 2n± and the peak height will 
decrease by 2p st — 1, where Pst = (l+n±)/(l + 2n±) is the 
ground state occupation number. These expressions are appli- 
cable at low temperatures, when the average thermal phonon 
number n t h = l/( ex P [^r/ksT] — 1) is small (n^ < 0.2 see, 
e.g., Ref!% For co r /27T = 8 GHz the Rabi peaks will be 
seen at 350 mK, but are negligible at IK (Fig. 2b). Actu- 
ally, at higher temperatures, when n^ > 1, the lowest dressed 
states become saturated and in addition to these there will ap- 
pear two broadened peaks (inside the Rabi doublet, Fig. 2a), 
roughly at frequencies ±g[\/n + 1 — y/n\ , n w n^, so that at IK 



these peaks will dominate over the Rabi peakJ^. These peaks 
could be a first quantum signature of acceptor-resonator cou- 
pling, when experiments are performed at the relatively high 
temperature of IK. 

Measurement via photons and standard techniques. The 
phonon probing can be performed using a recently proposed 
phonon-to-photon translator (PPT) 30 that can be realized on 
the same patterned Si nanomembrane with a simultaneous 
photon/phonon bandgap structure, Fig. lc. The translator al- 
lows for optical measurement techniques to be applied to 
phononic components. We show in Fig. 2e a schematic of the 
experiment, designed to measure the phonon cavity field am- 
plitude (b) via a heterodyne optical measurement 28 . To be 
able to scan around the mechanical resonance co r , one needs 
optical frequency resolution (at CO/27T w 200 THz 30 ) better 
than the dressed state peak width, (r+ k) /2n w 30 - 150 kHz, 
(at co r /27l c± 4 — 8 GHz, Table 1). This can be achieved by 
stabilization of the external photon probe to a bandwidth 
of tens of kHz. In the on-chip implementation of the PPT 
interface no photon should be allowed to enter the phoni- 
ton system, to avoid ionization of the acceptor (Methods). 
However, the estimated ionization cross section is small, 
tfphot ~ 8.6 x 10 -23 m 2 , since interband transitions are forbid- 
den: h(0 ~ 0.82eV < indirect bandgap in Si, and due to the 
final state energy suppression, ~ l/Ej- 5 . (For 10 photons in 

the cavity volume dX 2 one gets an acceptor lifetime of 12/is, 
assuming maximum photon- acceptor overlap.) Note that ion- 
ization can be further suppressed in regions of low photon in- 
tensity. 

The strong dispersive coupling regime is also reachable. 
The dispersive coupling, % = g 2 /A ar , ranges from 0.04 MHz 
to 2. 14 MHz for the frequencies presented (Table 1). Since the 
conditions A ar > 10 g (dispersive regime) andg 2 /A ar r re i ax > 1 
(good resolution of phonon number peaks) can be well ful- 
filled in the GHz range, this would allow an experiment on 
resolving the number states in the phononic cavity. In the dis- 
persive regime one can apply two tones (as in circuit- QEE^H): 
here, one tone is a phonon probe at C0j, slightly detuned from 
the resonator frequency co r (Fig. 2e). The second (spectro- 
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FIG. 2. Experimental signatures of a phoniton system, (a) Rabi splitting in the strong resonant coupling regime as a function of a phonon 
probe signal sweep, (b) Rabi peaks vs. temperature (approximate). At > 1 broad peaks will appear inside the Rabi doublet (not shown) due 
to transitions to higher dressed states (see text), (c) "Anticrossing picture" at qubit-phonon cavity detuning; for large positive detuning (upmost 
curve), the left resonance is dispersively shifted by —g 2 /A ar while the other (detuned) resonance is strongly suppressed as sin@o (see text), 
(d) An alternative experiment to determine the Rabi splitting: the absorption spectrum as a function of qubit detuning (via a magnetic field 
sweep), (e) Schematics of the standard optical heterodyne reflection experiment of the phoniton system utilizing a high-efficiency coherent 
photon-to-phonon translator. 



scopic) tone, at co^, is driving some of the dressed transi- 
tions around the acceptor frequency co a via direct microwave 
(MW) electric driving, that is similar to a MW manipulation 
of individual nitrogen- vacancy centers 13 (for B:Si a stronger 
coupling to electric MWs is expected comp ared to NVs since 
the acceptor dipole coupling is pe — Thus, one 

could measure the reflection by the phononic bandgap res- 
onator of the first (phononic) tone as a function of the MW 
spectroscopic tone (0 sp , and eventually observe the fine spec- 
tral structure of the dressed cavity-acceptor system 31 . We note 
that other measurement approaches are possible via, e.g., hole 
transport 2032 , STM probe spectroscopy 33 , or direct single- 
phonon detection^. 

Discussion and applications. We have introduced a sys- 
tem that allows for on-chip manipulation of coherent acoustic 
phonons via coupling to acceptor qubit states in a nanome- 
chanical cavity. Hybridization of the phonon- acceptor sys- 
tem and strong dispersive coupling are both possible with 
comparable parameters to circuit-QEB 21 and far surpassing 
semiconductor quantum dot QED^. The phoniton compo- 
nent can be incorporated in more complex systems such as 
with phonon-photon interfaces to photonic£^QI, and in arrays 
of other phoniton systems for engineered many -body phonon 
device d 14 ! 35 ! From the perspective of qubits, the isolated ac- 
ceptor provides a potentially robust two-level system for quan- 



tum information processing. The phoniton component de- 
scribed here offers an avenue for phonon dispersive readout of 
acceptor qubit states and the potential for spin qubit-to-photon 
conversion in silicon - which would be a seminal accomplish- 
ment. 



I. METHODS 

Origin of strong acceptor-phonon coupling. The role of 
an atom can be taken by a single impurity in a host crys- 
tal, e.g. in S™. The impurity-acoustic phonon interactional, 

ph( r ) = ^ij^ij ^ij ( r )' ma y l eac * to a stron g coupling regime 
(g > r re i ax ,K) even for cavity effective volume of few ten^ 
of ~ X 3 , since the deformation potential matrix elements are 
large: 15 (\\f s f\Dij\\\f s ) ~ eV. Qualitatively, the large coupling 
can be traced from the much smaller bandgap (~ eV) in the 
"Si-vacuum" as compared to QED (~ 10 6 eV). 

Hole val ence bands in Si. Holes in Si require a richer phys- 
ical pictur e 15 ! 16 ! (compared to positrons in QED). The 4-fold 
degeneracy (Fig. la) at the top of the valence band (neglect- 
ing heavy-light hole splitting) corresponds to propagation of 
particles of spin J = 3/2, reflecting the Tg representation of 
cubic symmetry. Relatively large spin-orbit coupling implies 
a 2-fold degenerate band representation), split-off by an 
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energy gap A S o ^ 45 For shallow acceptor centers in 

Si (e.g., B, Al, In, etc) the ionization energy, Ea ~ Aso: thus, 
all valence bands will play a role in the acceptor states. Still, 
the lowest acceptor states remain 4-fold degenerate since the 
acceptor spherical (Coulomb) potential does not change the 
cubic symmetry of the host crystal. 

Engineering the qubit levels. The role of electric field. The 
degeneracy of the ground state can be lifted by external mag- 
netic field via the Zeeman type interaction Hh (Fig. lb), and 
via mechanical strain, Eq.Q. The related Hamiltonians are 
invariants of the cubic symmetry group Oh = x I and time 
reversa l 15 ^ 16 ! and are constructed from the momentum opera- 
tor, k a 



IJL . 

i dx a 



e -A a (or the corresponding fields 15 ) and the 
spin-3/2 operators / a , a = x, y,z. 

For a relatively weak electric field E the linear Stark effect 
is possible: 



Pe_ 
V3 



{E X {jyJ z } + +Ey{j z J X }+ +E z {j X Jy} + ) , (4) 



since an ion impurity actually reduces the cubic symme- 
try (Td x I) to Td (and thus, there is no invariance under 
inversion) 15 . The ground state splits to two doubly degen- 
erate levels; however, the He does not commute with J z for 
any direction of the field E, leading to mixing of the Zeeman 
states. The latter can be useful to switch on/off the phonon 
coupling of the alternate qubit, { | — 1/2) , 1 1/2) }, provided the 
splitting 6Ee = 2pe\E\ is of the order of that due to stress, 
e.g., in the GHz range. The transition electric dipole moment, 
Pe, can be extracted from experiments: bulk dielectric ab- 
sorption measurements 1 ^ give pe — 0.26 D, with D = 3.336 x 
10 -30 C m being the Debye unit for e.d.m. (this is supported 
by single acceptor transport experiment d^3). Thus, a splitting 
of 1GHz requires an electric filed |£|iGhz — 3.85 x 10 5 V/m, 
achievable in nanodevicei^l Note, however, that increasing 
the field (splitting) exponentially decreases the qubit life time 
due to acceptor ionization: for $E = 1 GHz the life time is 
lion « 12 s, while for ~ 1.26 GHz it is Xi on « 12ms, etcPSl. 

These numbers show that there is an experimental "win- 
dow" for the alternate qubit, {| — l/2),|l/2)}, introduced in 
the main text. For example, for a qubit (Zeeman) splitting 
of dEn = 1 GHz and a strain splitting 8£ £ = 1 .43 GHz (ratio 
of ru = =0.7) the coupling factor reaches the maximal 
value of f(r h , r e ) ~ 0.4 for 8E E = 1 GHz, i.e., r h = 0.7. Anal- 
ogously, for a qubit splitting of 8£# = 2 GHz this electric field 
splitting leads to the same coupling factor of 0.4, giving a pos- 
sibility for a strong acceptor-phonon coupling, Eq.Q, and a 
relatively long (static field) ionization life time. 

Strong coupling of acceptor to confined acoustic phonon. 
We account for the acceptor coupling to a quantized phonon 
field starting from the Bir-Pikus Hamiltonian, derived for a 
uniform classical strain field, Eq.Q. For low-energy acous- 
tic phonons the interaction Hamiltonian, H^, has the form 



ex- 



of Eq.dl ), with the strain operator 8^(r) = \ -r ^ 
pressea via the quantized mechanical displacement: n(r) = 
(uqoir) bqo + u qa( r ) bqa^ ■ The mode normalization is 
f d 3 ru qa (r)u qa (r) = 2p ^ , so that b\ Q creates a phonon in 



the mode (q, a) with energy h(ti qa (p is the material mass den- 
sity) in a mode volume V. The vector q denotes a collective 
index of the discrete phonon mode defined via the phonon 
cavity boundary conditions and mode volume V . Simi- 
lar to cavity QEB 28 , the phonon- acceptor coupling hg s qQ = 
(s f ; qa | H p h | s) enters in a Jaynes-Cummings Hamiltonian (see, 

e.g. RefPS): H g « hg s q ^ (o+, s b q p + cT^<U) ' where we onl y 
retain the resonant cavity phonon, and G^, s = \s')(s\ refers to 
the relevant acceptor transition. 

Phonon-photon translator. Such a device is based on op- 
tomechanical non-linearities that couple in the same bandgap 
cavity (see Fig. 2e) two photon modes (a,a p ) and a phonon 
mode b, via optomechanical coupling h om 30 . For photons 
in the near-infrared range (k opt w 1500nm) the PPT allows 
one to couple a quantum optical input/output channel (of fre- 
quency CO/27T ~ 200 THz) to a phonon channel (with COj/27T ~ 
4 — 8 GHz), and the coupling between the fields is enhanced by 
the auxiliary photon pump channel, pumping at the sideband 
resolved frequency co^ = CO — co m — A (at pump detuning A = 
it is at resonance with the red sideband of mode CO). The co- 
herent nature of the photon-to-phonon translator is described 
by the effective beam-splitter type Hamiltonian 36 : 



A^S + Go 



(d^b + dP^J 



(5) 



where G om ^ h om Eo is the enhanced effective coupling, pro- 
portional to the pump field amplitude, Eq. The weak coupling 
regime, G om < K° pt is needed to avoid total optical reflec- 
tion, and optimal translation (close to 100%) takes place at 
a matching condition 30 G om y/^nech^opt ^mech^ K o P t are the 
couplings of the PPT to respective waveguides). 

Acceptor ionization via optical photons. Since the photon- 
to-phonon translator is realized on the same Si nanomembrane 
(implying a simultaneous photonic/phononic bandgap struc- 
ture) it is natural to ask how 200 THz photons may affect the 
qubit lifetime when they reach the acceptor. The correspond- 
ing photon energy of 0.82 eV is less than the indirect bandgap 
in Si (A£ gap = 1.1 eV) and thus interband transitions are not 
possible. Thus, one considers an "ionization process" of a 
bound hole going to the continuous spectrum, an analog of 
the ionization of an (anti)hydrogen atom (the corresponding 
cross section is thus re- scaled). Correspondingly, one uses 
the re-scaled values: a free hole mass yua — 0.23 m e in Si, an 
effective Bohr radius af f = 2 [4 7r e f ]e a ' w * m me acce P tor i° n_ 
ization energy for B:Si, Ea ~ 0.044 eV, and screening factor 
Z ~ 1.4. The total cross section is: 



Gphot 



3271 



h 6 



1 



3 cV2n^m 3 A [af] 5 E 2 f 5 (E f + E A ) [4TOoe r ] ' 



(6) 



where Ef = to — Ea is the final (free) hole energy, and c = 
Co I V^r is the speed of light in Si (zf 1 ~ 1 1 .9). Since Ea^Ej, 
the total cross section is suppressed as oc l/£" 3 - 5 . Given n c 
photons in cavity volume V the acceptor life time is x p h t = 
2V I (ft c cc>phot) for a maximum photon-acceptor overlap. 
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